An approach of dynamic mesh adaptation for simulating 3‐dimensional unsteady moving‐immersed‐boundary flows

In this paper, we present an approach of dynamic mesh adaptation for simulating complex 3‐dimensional incompressible moving‐boundary flows by immersed boundary methods. Tetrahedral meshes are adapted by a hierarchical refining/coarsening algorithm. Regular refinement is accomplished by dividing 1 tetrahedron into 8 subcells, and irregular refinement is only for eliminating the hanging points. Merging the 8 subcells obtained by regular refinement, the mesh is coarsened. With hierarchical refining/coarsening, mesh adaptivity can be achieved by adjusting the mesh only 1 time for each adaptation period. The level difference between 2 neighboring cells never exceeds 1, and the geometrical quality of mesh does not degrade as the level of adaptive mesh increases. A predictor‐corrector scheme is introduced to eliminate the phase lag between adapted mesh and unsteady solution. The error caused by each solution transferring from the old mesh to the new adapted one is small because most of the nodes on the 2 meshes are coincident. An immersed boundary method named local domain‐free discretization is employed to solve the flow equations. Several numerical experiments have been conducted for 3‐dimensional incompressible moving‐boundary flows. By using the present approach, the number of mesh nodes is reduced greatly while the accuracy of solution can be preserved.     

Experimental and numerical study of an aquaculture net cage with floater in waves and current

The mooring loads on an aquaculture net cage in current and waves are investigated by dedicated model tests and numerical simulations. The main purpose is to investigate which physical effects are dominant for mooring loads, and in this respect, to investigate the validity of different rational hydrodynamic load models. Also structural and numerical aspects are investigated. The model tests are performed to provide benchmark data, while the numerical model is used to study the effect and sensitivity of different load models and parameters.Compared to a realistic aquaculture plant, the total system is simplified to reduce the complexity. The system does, however, include all the four main components of an aquaculture plant: net cage, floater, sinker weights and moorings. The net cage is bottomless, flexible and circular. It is attached to a circular, elastic floater at the top and has 16 sinker weights at the bottom. The system is nearly linearly moored with four crow feet mooring lines.The loads are measured in the four mooring lines. A systematic variation of current only, wave only as well as combined current and wave conditions is carried out. The numerical simulation results are first benchmarked towards the experimental data. The mean loads in general dominate over the dynamic part of the loads in combined current and waves, and they significantly increase in long and steep waves, relative to current only. Next, a sensitivity study is carried out. A rigid floater significantly alters the loads in the mooring lines compared to a realistic, elastic floater. The theoretical model for the wave matters. The mooring loads are rather insensitive to a majority of the parameters and models, in particular: frequency dependent added mass of the floater and nonlinear restoring loads. It seems not to be necessary to represent the net cage with a very fine numerical mesh.     

船舶与海洋平台结构冰载荷的高性能扩展多面体离散元方法

船舶与海洋平台结构的冰载荷是寒区海洋工程结构物设计中的关键参数,而离散元方法是有效计算结构冰载荷的重要手段. 本文采用基于闵可夫斯基和原理的扩展多面体离散元方法模拟船舶与海洋平台结构的相互作用过程. 其中,构造扩展多面体的近似包络函数并建立了基于优化模型的快速接触搜索算法;考虑单元间粘结作用的刚度软化过程建标识码元间的粘结-破碎模型. 同时,发展了 CPU-GPU 协同异构环境下的高性能并行算法. 为分析海冰与海洋结构作用中的冰载荷,采用ISO标准验证了扩展多面体离散元分析结构冰载荷的准确性. 采用离散元方法计算了船舶结构的冰载荷,研究了船舶结构表明的线载荷分布特点,并采用船舶结构冰阻力经验公式验证了计算结果的合理性. 采用离散元方法计算了平整冰区与多桩腿平台结构的相互作用,分析各桩腿上的冰载荷特点. 针对碎冰区的海冰管理过程,采用离散元方法分析了船舶结构绕行过程中的船舶和海洋平台结构冰载荷. 本文方法可有效应用于海洋结构冰载荷分析,能为极地船舶与海洋平台结构的设计和安全运行提供科学的分析手段.      

基于Voronoi切割算法的碎冰区构造及离散元分析

离散元方法广泛应用于海冰，特别是碎冰区的动力过程及其对海洋结构作用过程的数值模拟。为构造碎冰区中的冰块几何特性，基于二维Voronoi图方法对计算域进行随机切割以生成碎冰区中冰块的几何形态，并采用球体单元对每个碎冰块单元进行填充，从而确定碎冰区的初始分布场。在采用Voronoi图进行碎冰区构造时，可对冰块尺寸、几何形态和密集度等海冰参数进行设定。为确定冰块的不同几何规则度，综合采用排斥法和扰动法以定量地控制碎冰块几何形态从完全随机分布到规则分布的连续变换。为分析不同几何规则度下碎冰块的几何特性概率分布规律，对计算域内冰块的面积和边数等参数进行统计分析，从而可更合理地参数化控制初始冰场中碎冰块的几何特性。在此基础上，本文基于粘接-破碎的球体离散元方法对不同冰况下锥体结构的冰荷载进行了数值计算，讨论分析了碎冰区的海冰密集度、冰块面积和几何规则度对冰载荷的影响。     

水面舰艇编队的最佳防空队形探讨

最佳防御队形以编队对来袭导弹的可探测面积尽可能大为前提,并以"抗饱和攻击能力"为衡量标准.各方向的可拦截批次受两方面因素限制:一是来袭导弹被发现时其与指挥舰的距离,一般距离越大,防御准备就越充分,可拦截批次就越大;二是护卫舰到来袭导弹轨迹的距离,一般距离越小,单次拦截时间就越短,可拦截批次就越大.定义以概率1可拦截批次最小的方向为最危险方向,经计算初始队形各方向可拦截的批次不等,通过"削峰补谷"的方式予以均衡和优化.若以拦截批次的期望为标准,最危险方向与以概率1可拦截的批次为标准的结果相同.如果得到空中预警机的信息支援,在最危险方向上编队就可更早地对距离指挥舰148.4km远的导弹发起拦截,增大编队的抗饱和攻击能力,但由于防空导弹射程限制,预警机提供的信,息支援无法得到充分利用,此时限制编队抗饱和攻击能力的主要矛盾转向防空导弹的射程.     

Allostatic load as a complex clinical construct: A case‐based computational modeling approach

Allostatic load (AL) is a complex clinical construct, providing a unique window into the cumulative impact of stress. However, due to its inherent complexity, AL presents two major measurement challenges to conventional statistical modeling (the field's dominant methodology): it is comprised of a complex causal network of bioallostatic systems, represented by an even larger set of dynamic biomarkers; and, it is situated within a web of antecedent socioecological systems, linking AL to differences in health outcomes and disparities. To address these challenges, we employed case‐based computational modeling (CBM), which allowed us to make four advances: (1) we developed a multisystem, 7‐factor (20 biomarker) model of AL's network of allostatic systems; (2) used it to create a catalog of nine different clinical AL profiles (causal pathways); (3) linked each clinical profile to a typology of 23 health outcomes; and (4) explored our results (post hoc) as a function of gender, a key socioecological factor. In terms of highlights, (a) the Healthy clinical profile had few health risks; (b) the pro‐inflammatory profile linked to high blood pressure and diabetes; (c) Low Stress Hormones linked to heart disease, TIA/Stroke, diabetes, and circulation problems; and (d) high stress hormones linked to heart disease and high blood pressure. Post hoc analyses also found that males were overrepresented on the High Blood Pressure (61.2%), Metabolic Syndrome (63.2%), High Stress Hormones (66.4%), and High Blood Sugar (57.1%); while females were overrepresented on the Healthy (81.9%), Low Stress Hormones (66.3%), and Low Stress Antagonists (stress buffers) (95.4%) profiles. © 2015 Wiley Periodicals, Inc. Complexity 21: 291–306, 2016     

Radial sample load distribution under overload conditions: Analytical scale columns

The homogeneity of the sample load across the radial cross section of analytical scale columns was determined when operating under overload conditions. The study was performed using active flow technology columns operating in parallel segmentation mode. The outlet segmentation ratio was varied to enable different volume fractions of mobile phase, and thus sample, to elute from the peripheral and central flow regions of the column. The amount of solute exiting the peripheral and radial central exit ports was determined as a function of the flow segmentation ratio. The experimental data using an analytical scale column with dimensions, 100?×?4.6?mm, indicated that the sample load distribution was essentially uniform as a function of the column radial cross section.     

An ALE finite element method for a coupled Stefan problem and Navier–Stokes equations with free capillary surface

In this article, an ALE finite element method to simulate the partial melting of a workpiece of metal is presented. The model includes the heat transport in both the solid and liquid part, fluid flow in the liquid phase by the Navier–Stokes equations, tracking of the melt interface solid/liquid by the Stefan condition, treatment of the capillary boundary accounting for surface tension effects and a radiative boundary condition. We show that an accurate treatment of the moving boundaries is crucial to resolve their respective influences on the flow field and thus on the overall energy transport correctly. This is achieved by a mesh‐moving method, which explicitly tracks the phase boundary and makes it possible to use a sharp interface model without singularities in the boundary conditions at the triple junction. A numerical example describing the welding of a thin‐steel wire end by a laser, where all aforementioned effects have to be taken into account, proves the effectiveness of the approach.Copyright © 2012 John Wiley & Sons, Ltd.     

The Calculus of Moving Surfaces and Laplace Eigenvalues on an Ellipse with Low Eccentricity

Our goal is to demonstrate the utility of the calculus of moving surfaces (CMS) in boundary variation problems. We discuss the relative advantages of the CMS compared to the alternative approach of interior variations. We illustrate the technique by calculating the two leading terms of a power series for the Laplace eigenvalues on an ellipse with semi-axes 1 + a and 1 + b, where a and b are small. We compare the CMS estimates with those obtained by the conventional finite element method with Richardson extrapolation. The comparison confirms the cubic rate of convergence for the CMS estimates.     

Moving grid method for numerical simulation of stratified flows

We develop a second‐order accurate Navier–Stokes solver based on r‐adaptivity of the underlying numerical discretization. The motion of the mesh is based on the fluid velocity field; however, certain adjustments to the Lagrangian velocities are introduced to maintain quality of the mesh. The adjustments are based on the variational approach of energy minimization to redistribute grid points closer to the areas of rapid solution variation. To quantify the numerical diffusion inherent to each method, we monitor changes in the background potential energy, computation of which is based on the density field. We demonstrate on a standing interfacial gravity wave simulation how using our method of grid evolution decreases the rate of increase of the background potential energy compared with using the same advection scheme on the stationary grid. To further highlight the benefit of the proposed moving grid method, we apply it to the nonhydrostatic lock‐exchange flow where the evolution of the interface is more complex than in the standing wave test case. Naive grid evolution based on the fluid velocities in the lock‐exchange flow leads to grid tangling as Kelvin–Helmholtz billows develop at the interface. This is remedied by grid refinement using the variational approach. Copyright © 2012 John Wiley & Sons, Ltd.